The Radio-FIR correlation: Is MHD Turbulence the 

Cause? 
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Abstract 

The radio - far infrared correlation is one of the tightest correlations found 
in astronomy. Many of the models explaining this correlation rely on the associ- 
ation of of global magnetic field strength with gas density. In this letter we put 
forward that the physical reason for this association lies within the processes of 
magnetohydrodynamic turbulence. 
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1 Introduction 

One of the more extraordinary correlations in astronomy is that between the far in- 
frared (FIR) and radio continuum of galaxies. A correlation betw een the IR emis- 
sion and radio was first noted in Seyferts by Ivan der Kruitl (|l97lh . it was not un- 
til the results of th e IRAS mission were analysed that universality of this correla - 
tion was discovered ( Dickev fc Salpeter . 1984 : de Jong et al. . 19851: Helou et al. . 1985| ). 
This linear correlation spans ~ 5 orders of magnitude with less than 50% disper- 



sion (|Yun. Reddv. fc Condon . 2001: IWunderlich. Wielebinski. fc Kleinl . 119 87). making 
it one of the tightest correlations known in astronomy. It appears to be followed b y 
all galaxies with ongoing star-formation and without a dominant AGN 
The correlation not only holds on global scales but is als o found to hold within the disks 
of galaxies, down to scales on the order of a few 100 pc ( Beck fc Golla . 1988| ). What is 
most extraordinary about this relationship is that it couples a purely thermal process 
in the IR (the re-emission of UV radiation by dust grains) with a non-thermal process 
at radio wavelengths (the synchrotron radiation of relativistic electrons). Typically, at 
1.4 GHz, the non-ther mal emission dominates by at least an order of magnitude over 
the free-free emission (|CondonL Il992t iNiklaj . Il997t iDooita et all 120021 ) , and hence the 
contribution of free-free emission to the correlation is minimal. 
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T he radio - FIR c orrelation is not a simple mass-scaling (richness) effect ((Wunderlich fc Klein , 

1988| Xu et al. . 1994 ). and hence must be explained by some form of direct coupling be- 



tween the IR emitting dust and the synchrotron emission o f cosmic rays. Star formatio n 
is generally accepted to be responsible for the correlation ( Wunderlich fc Klein . 1988| ). 
The basic scenario is as follows: As the massive, hot stars are born in star-forming 
regions and live their short lives, the UV radiation they emit heats the surrounding 
dust, which then re-emits in the IR. This radiation also provides the energy for the 
thermal (free-free) radio emission. When these hot stars reach the end of their brief 
existence, the resulting supernova shocks create energetic cosmic rays which are be- 
lieved to be responsible for the non-thermal radio emission. However, as it stands this 
model has too many steps and too many param eters to explain the tightness of the 
correlation (see the illustration bv lEkersl (jl99lh ). Thus more complex models have 
been put forward, although there is still no consensus on the cause of the radio - FIR 
correlation, or an adequate explanation for its tightness. In several of these models a 
large part of the explanation of the correlation was related to the association of global 
magnetic field strength with gas density through an energy equipartition scheme. In 
this letter we put forward a physical mechanism for this assumption adopted in these 
previous models. We show here that the cause for this relationship may lie within the 
processes of magnetohydrodynamic (MHD) turbulence. 

In section 2 we discuss the previous theoretical models that have been previously 
proposed to explain the correlation. In section 3 we discuss MHD turbulence and how 
it naturally provides the required relationship between magnetic field strength and gas 
density. In section 3 we discuss how this relation provides a possible basis for the 
correlation, and how it fits into the previous models, with the concluding remarks in 
section 4. 



2 Previous Theoretical Models 



As mentioned in the introduction, recent star-formation provides a basis for most theo- 
retical models of the FIR-radio corre l ation. One of the earliest theories, the 'optically- 
thick' or 'calorimeter' theory (|Volkl . Il98flt IVolk fc XiJ . 119941 iLisenfeld Volk. fc Xu . 
1996), assumed three things. First, that all the FUV 1 radiation from massive stars is 
absorbed by the dust grains within a galaxy. Second, that the energetic cosmic rays 
produced by the supernova explosions of these stars lose most of their energy within 
the galaxy due to synchrotron and inverse Compton process. As both these processes 
are proportional to the number of massive stars, these calorimetric assumptions lead 
to the correlation. Finally, the tightness of the correlation is provided by the third 
assumption that the energy density of the interstellar radiation field, U ra d, is in a con- 
stant ratio with the magnetic field energy density, Ub- In other models which do use 
the first assumption ('optically thin' models), including ours, there is the generally im- 
plicit assumption that there is a direct linear relation between the gas de nsity and dust 
density, an assumption which is nonetheless supported by observations (|Xu &: Helou . 



1 FUV radiation is that which lies between « 5 eV and w 13.6 eV 
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1996). 



An alternative theory put forward by Helou fc Bicav ( 1993| ) assumes the opposite 
extreme, an 'optically-thin' model, in which the cosmic rays and UV photons both have 
high escape probabilities. To provide the correlation in these optically and CR thin 
galaxies they have two assumptions. Firstly that the UV or 'dust-heating' photons and 
the radio emitting cosmic rays are created in constant proportion to each other, which 
is again related to star formation. Secondly, the tightness of the correlation is provided 
by a local coupling between the magnetic field strengt h and the gas density. 

A challenge to both these theories is the model bv lNiklas fc Beckl ()1997l ). In this 
work they argue that observations indicate that within most galaxies, CR electrons 
lose very little energy before they escape. These same galaxies are optically thick to 
UV photons, thus both the calorimetric and optically-thin models are not supported 
by the observations. In the Niklas & Beck model, the controlling factor they put 
forward for the correlation is the volume density of the gas. They assume that both 
the star formation rate (and thus dust heating) and the magnetic field strength (which 
determines the synchrotron emission) depend upon the gas volume density and hence 
the correlation. 

Of course there have been other approaches, such as that of iBettens et al. (1993) 
who looked at cosmic ray driven chemistry in molecular clouds to explain the correla- 
tion. However, in their case they did not take account of the strength of the interstellar 
magnetic field, an important factor in the correlation. 



A different way to look at the correlation was suggested bv lHoernes. Berkhuiisen. fc Xu 
(1998) when examining M31. They decomposed both the radio and the FIR into two 
parts. The radio was decomposed into a thermal radio component and a non-thermal 
component, while the FIR was decomposed into a warm component, associated with 
Hn regions, and a cool component, associated with the diffuse (cirrus) clouds. They 
found both a correlation between the thermal radio a nd warm FIR and a co r relatio n 
between the non-thermal radio and the cool FIR (like IXu. Lisenfeld. fc Volk ( 1994 )). 
The first correlation is easily understandable as both components are associated with 
Hn regions. The second correlation is not so easy to understand, as the cool FIR 
component of M31 is not believed to be predominantly heated by massive stars. The 
local low mass stars provide most of the heating, yet these are not the progenitors 
of supernovae which provide the CR electrons necessary for the non-thermal radio. 
Thus in M31 at least, the correlation of non-thermal radio - cool FIR must be due to 
factors other than star formation. What Hoernes. Berkhuiisen. fc Xu ( 1998| ) suggest is 
that the correlation is due to the close coupling of the magnetic field strength and gas 
density. Assuming that the gas density is proportional to the dust density then the 
synchrotron emission (via the magnetic field) is directly related to the FIR emission 
(due to the local dust density). 

It is clear that most of these theories require some form of coupling between the 
magnetic field and gas density. It is usually asserted that 'equipartition' provides the 
coupling, yet the way in which equipartition can arise is left open. In the next section 
we explore whether MHD turbulence can provide the necessary coupling mechanism. 
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3 Turbulent Coupling 



3.1 MHD Turbulence 

Turbulent motions are observed in most astrophysical fluids and since magnetic fields 
are undoubtedly present in such fluids, magn etohydrodynamic (MHD) turbule nce is 
an important field of study in astronomy (see Cho. Lazarian. fc Vishniac (20,02) for a 
review of MHD turbulence). 

In most astrophysical plasmas, the magnetic Reynolds number (R m = L5V/ rj, where 
L is the size of the system, 5V is the the r.m.s. velocity, and rj is the magnetic 
diffusion) easily exceeds 10 10 and the usual expectation is that magnetic field is frozen 
into the gas in such systems. The velocity field advects and stretches magnetic field 
lines and magnetic field exerts pressure and tension forces on velocity fields. With 
these additional considerations, MHD turbulence is generally different from the pure 
hydrodynamic case. 

On the scales over which the radio - FIR correlation holds, it is reasonable to 
assume that the over-all magnetic field is weak. In this weak/zero mean field regime, 
there are two main mechanism s for the generation of m agnetic field: the dynamo effect 
( Parkerl . E97i iMoffat A Il97ah and field line stretching (|Batchelorl . Il95(ih . The dynamo 



effect can amplify the mean magnetic field (-Bo); while field line stretching is responsible 
for amplification of the local mag netic field. 



Using MHD turbulent models, ICho fc Vishniac (2000a) have shown that even with 



out the dynamo effect, field line stretching can amplify the magnetic fields up to the 
level of energy equipartition. According to this model, the rate of field line stretching 
at the scale of energy injection or the largest energy containing eddies, L, is ~ SV/L, 
while the rate of magnetic back-reaction at the scale is ~ (SB/^AiTp)/L, where SV is 
the r.m.s. velocity, SB is the r.m.s. magnetic field strength, and p is the average density. 

Therefore, when SB/^4np <C SV, stretching is more effective than back- reaction, 
resulting in the growth of the r.m.s. field strength. The stretching is balanced by the 
back-reaction only when the local Alfven velocity (5B / y/Anp) is comparable with the 
local fluid velocity {SV). In other words, magnetic field cannot further grow by field 
stretching when the energy equipartition condition, p(SV) 2 ~ (SB) 2 /(4ir), has been 
reached. That is, when MHD turbulence reaches the stationary state, the magnetic 
energy density matches the energy density of the gas, SB / \fAixp ~ SV . In the weak 
field case the turbulent local magnetic field {SB) is larger than the external, mean 
magnetic field (B ). 

In principle, the dynamo can amplify the large-scale magnetic field (-Bo)- I n mean 
field dynamo theory (see Moffatt 1978; Parker 1979), turbulent motions at small scales 
are biased to create a non-zero electromotive force along the direction of the large-scale 
magnetic field. This effect (the 'a-effect') works to amplify and maintain large-scale 
magnetic fields. Whether or not this effect actually works depends on the structure 
of the MHD turbulenc e , espe cially on the mobility of the field lines. For example, 
Vainshtein fc Cattaneol (jl992| ) have argued that when equipartition between magnetic 



and kinetic energy densities occurs at any scale larger than the dissipation scale, the 
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mobility of the field lines and the a effect will be greatly reduc ed. However the na- 



ture and degree of this suppression is a controversial issue (see Gruinov fc Diamondl 



( 19941) :ICattaneo fc Hughel (Il996l) Packman fc Field! (|2000MVishniac fc Choi f)200lh : 



see Vishniac. Lazarian fc Cha t 20021 ) for a review). Nevertheless, when the mean field 
grows to the value similar to energy equipartition (£> / y/knp ~ SV), mobility of mag- 
netic field lines is greatly reduced and ther e fore th e dynamo can no longer operate. Nu- 
merical simulations (e.g. ICho fc Vishniac ( 2000bl )) show that, when Bo/^Anp ~ 5V, 
there also exists (almost exact) energy equipartion between random magnetic and 
turbulent kinetic energy. 
5B/^ Wp ~ 5V. 

(B ~ a/-Bq + (SB) 2 ), as well as the fluctuating one (5B), stays at the equipartition 
value: B/y/Arfp ~ \pl 5V. 

Thus, regardless of the strength of the external (galactic) magnetic field, the local 
magnetic field will stay in equipartition with the gas; 



From this observation, we can assume that B / 'y/Anp ~ 
Therefore, when the dynamo operates, the local magnetic field 



5V ~ B/^fp. 



This provides the coupling between the gas density and the magnetic field which is a 
necessary, but not sufficient, condition for the operation of the radio - FIR correlation. 



3.2 Magnetic field and Gas Density Coupling 



When the velocity dispersion of quiescent gas within the disk of galaxies is measured it is 
found to be remarkably constant. HI observations o f the LMC found the velocity disper- 
sion l ies within 6.8 and 7.7 km s _1 across the galaxy ( Kim. Dopita. Stavelev-Smith. &: Besselll . 
Il999l ). Similar observations of the nearly face-on galaxies NGC 628 and NGC 3938 also 
show approximately constant velocity dispersion ac ross the face of the gal a xies in both 
CO (NGC 628: 6 kms" 1 , NGC 3938: 8.5 kms 1 ICombes fc Becauaertl l|l997h ) and 
HI (NGC 628: ~ 9 kms" 1 (IShostak fc van der KruitL Il984h . NGC 3938: ~ 10 km s" 1 
( van der Kruit fc Shostak . 19821 ) ). This trend of a velocity dispersion of ~ 10 kms -1 



continues beyond these three objects, with the dispersion being app reciably uniform be- 
tween many galaxi e s, inc luding both spirals and irregulars (see e.g. Sellwood k, Balbusl 
f)l999l ): lKamphuisl (|l993h . Though the reason behind this unusual constancy is not 
yet understood, we can still apply this fact to the situation of energy equipartition. 
Equation^ with a constant velocity dispersion then implies that: 



(2) 



B oc yp. 

This relation is consistent with other numerical MHD simulations. These sim- 
ulations of turbulent interstell ar gas show B oc Tf 1 with m found to be m ~0.4-0.6 
by several independent models (Padoan &: Nordlundl 19991: Kim. Balsara. &: Mac Low . 
2001t lOstriker. Stone, fc Gammid . 120011 ) . The relation between B and p shows that 
there is a tendency toward energy equipartition even within individual eddies, which 
may represent individual clouds. Note that the energy equipartition described in fcj.ll 
is for the entire system, which may represent units as large as an entire galaxy. Of 
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course, for this to represent an entire galaxy we must rely on the implicit assumption 
that the magnetic field within the halo is closely coupled to that within the disc. 

A similar correlation between B and p is also seen observationally within galaxies, 
with iBerkhuiisen fc et al.l ()1997l ) showing B oc p g m in M31 and the Milky Way with 
m ~0.3-0.7. Similarly, th ere lies an observation al correlation between B and p on 
global galactic scales, with iNiklas &: Beck! (|1997l ) using 43 galaxies to find a relation 
with a slope m = 0.48 ± 0.05. Note that in both cases equipartition between cosmic 
rays and magnetic field energy densities was assumed. 

Thus the relation betw een B and p is seen in observations and is also inferred in 
our numerical simulations ( Cho &: Vishniac , 2000a f]) , with MHD turbulence providing 
the mechanism for this coupling. 



4 The Connection to the Radio — FIR Correlation 

As discussed in §2, most models that provide an explanation for the radio - FIR 
correlation assume a coupling of the magnetic field and gas density of the form of eqn. 
El The physical mechanism put forward in the previous section justifies this assumption 
and thus provides further credence to these models. 

What the relation between B and p also provides is a direct coupling between two 
parameters which determine the amount of flux in the FIR and radio and hence relates 
directly to the correlation. If we igno re the effect of the the r mal radio emission - whic h 
is usually minimal on galactic scales (|Condonl . Il992t iNiklaj . Il997l : iDooita et all l2002h , 



then the remaining parameters of the correlation are the number density of the cosmic 
rays and the dust-heating UV radiation field. 

The density of cosmic rays is presumably determined by their production via Fermi 
acceleration processes occurring in strong shocks, mostly generated by stellar winds 
and by supernova explosions, their radiative lifetime, and the volume of the galaxy in 
which they can move. If, in a disk galaxy, £sf is the birthrate of young stars per unit 
area of the disk, and the characteristic scale-height of the cosmic rays is z cr , determined 
by the large-scale magnetic field configuration, then the number density of comic rays 
with be given by 

n CT = (p Esf / z CI . (3) 

where is a constant determined by the ratio of the characteristic acceleration efficiency 
and the radiative lifetime. These radiative losses will mostly occur in highly localized 
regions corresponding to the denser neutral clouds where the local magnetic field has 
been enhanced by MHD turbulence. The local synchrotron emissivity at frequency v 
is given by 

j v = f{a)kB^ 2 u-^ 2 (4) 

where the number density of the relativistic electrons with relativistic 7 has a power law 
distribution N( / y)~ a = Ary~ a , and B is the local magnetic field. The total density of the 
relativistic electrons is therefore n cr = kj^/ (a — 1) where a ~ 2.4 which corresponds 
to a frequency spectral index of synchrotron emission a = —0.7. Such a value for the 
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index is not unreasonable ( Stevens. Forbes, fc Norrisl 2002j ). though the value could be 



flatter or steeper than this. Therefore, the non-thermal emissivity scales as the product 
n CT B 17 , approximately, or: 

jv = <PB 1 - 7 SSF /Zcr (5) 

The local infrared emissivity scales as the local radiation field, which can be ap- 
proximated as: 

where if) and ( are scaling factors, is the characteristic scale height of the young stars 
and E* is the surface density of the cool, old stellar population. Note that the scale 
height of the old stellar population is larger than that of the young population but 
this factor is taken account within (. The local emissivity also depends on the local 
gas density and the radiation field- weighted mean opacity of the dust grains, (/«uv}, 
which though dependent upon the grain properties and radiation field, will vary only 
to a small extent. Thus, assuming that the dust emission is optically thin, we have 

Jfir = nip (kw) (Ssf +C£*)/z*. (6) 

Therefore, with our MHD turbulence, which provides the scaling relationship n oc B 2 , 
locally, we would expect that 

jv <f>Z* SSF . 0.15 



n^. (7) 

JFIR t/j (Kuv) 2cr(S S F +C^*) 

This provides the required radio-FIR correlation to the degree to which the product 
of the physical parameters on the right hand side of this equation remain constant. 
Of course this is only valid when the same volume element produces both the FIR 
and synchrotron emission. There are possibly other volume elements (such as in the 
halo) which might only produce radio emission. This means that equation[7|may not be 
relevant for the global FIR-radio correlation, only giving a quantitative analysis for the 
local correlation. Additionally, variations in these physical parameters with different 
galaxy types may also explain the non-linear ities and outliers in the correlation, such 



as starbursts (|Yun. Reddv. fc Condon! . 120011 1 . Further work needs to be done on the 



relation to the global correlation as well as the scatter in the parameters in equation 
13 and others, like the local synchrotron index, to determine their effects upon the 
correlation and whether further constraining mechanisms are needed to explain the 
tightness of the correlation. 

Thus, though we have explained one part of the radio - FIR correlation, more 
understanding is needed of the cosmic rays, radiation field heating the dust, and the 
processes discussed above before the exact reason for the correlation and its tightness 
can be fully comprehended. 

5 Conclusion 

Most models which try to explain the remarkable correlation between radio and FIR 
emission rely upon the association of magnetic field strength and gas density. MHD 



7 



simulations show that through the process of equipartition, the magnetic field and gas 
densities are coupled, with 



a relationship that is also indicated by observations within our own and other galaxies. 
This relationship provides a basis for the radio - FIR correlation by directly connecting 
one of the parameters of synchrotron emission with a parameter of FIR emission. While 
not explaining the relationship fully, it does put us a step closer to fully understanding 
this remarkable correlation. 
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